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ABSTRACT 

NGC 1360 is a large planetary nebula (PN) without an obvious shell morphology. We have analyzed Ha 
images and high-dispersion echelle spectra of NGC 1360 in order to construct spatio-kinematic models 
and to determine its density distribution. The best-fit model indicates that NGC 1360 is a prolate 
ellipsoidal shell whose major axis is twice as long as its minor axis and is tilted by 60° with respect to 
the line of sight. The large kinematic age of the shell, ^10,000 yr, and the low density of the nebula, 
<130 H-atom cm -3 , imply that NGC 1360 is an evolved PN and has begun to merge with the interstellar 
medium. The observed morphology and surface brightness profiles of NGC 1360 can be described well 
as a thick shell with a Gaussian radial density profile without a sharp inner edge, indicating a lack of 
on-going compression by a fast stellar wind. The FLIERs observed in NGC 1360 near the end of its 
major axis expand faster and are younger than the nebular shell. 

Subject headings: ISM: kinematics and dynamics - ISM: planetary nebulae: general - planetary 
nebulae: individual (NGC 1360) 



successfully reproduced by interacting-stellar-wind mod- 
els. 

There exist PNe that do not possess morphological fea- 
tures to suggest the presence of wind-wind interactions. 
NGC 1360 is such a PN. Its surface brightness does not dip 
deeply at the center or rise steeply at the limb to indicate 
a hollow shell structure. The lack of shell morphology in 
NGC 1360 makes it particularly interesting. Therefore, we 
have chosen this PN for a spatio-kinematic analysis in or- 
der to determine its physical structure and to investigate 
its evolutionary status. This paper reports our observa- 
tions and analysis of NGC 1360. 

2. OBSERVATIONS 

2.1. Direct Images 

Narrow-band Ha images of NGC 1360 were obtained on 
1999 October 8-9 with the 1 m telescope at Mount Laguna 
Observatory (MLO). The central wavelength and FWHM 
of the Ha filter are 6563 A and 17 A, respectively. The 
detector was a Loral 2K CCD with a pixel size of 15 /zm, 
corresponding to an image scale of 0'.'4 pixel -1 and provid- 
ing a field of view of 13'3 x 13f3. The angular resolution 
was ~ 2'.'8, as determined from the FWHM of field stars 
in the images. A total of 10 images, each with an exposure 
time of 600 s, were combined to form the final image. This 
image is shown in Figure 1. 

We have also retrieved archival Hubble Space Telescope 
images of NGC 1360 taken with the Wide Field and Plane- 
tary Camera 2 (WFPC2) and the F656N (Ha) and F658N 
([N II]) filters. These observations were obtained on 1994 
February 15, during the "pre-cooldown" period when the 
WFPC2 was operating at a temperature of —77 °C (rather 

Current address: Instituto de AstrofTsica de Andaluci'a, Consejo Superior de Investigaciones Cienti'ficas (CSIC), Apartado Correos 3004, 
E- 18080, Granada, Spain. 

2 Visiting Astronomer, Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatories, operated by the Association of 
Universities for Research in Astronomy, Inc. (AURA) under a cooperative agreement with the National Science Foundation. 

1 



1. INTRODUCTION 

Planetary nebulae (PNe) exhibit a wide range of mor- 
phologies. While the overall shape of a PN may be round, 
elliptical, or bipolar, most are limb-brightened, indicating 
a lower density in the interior, or the existence of a cen- 
tral cavity. It has been suggested that PNe are formed 
by dynamical interactions between the current fast stellar 
wind and previous slow wind (Kwok, Purton, & Fitzgerald 
1978). Hydrodynamic calculations of interacting-stellar- 
wind models have indeed been successful in reproducing 
the hollow-shell structure and overall shapes of PNe by 
varying the pole-to-equator density gradient in the Asymp- 
totic Giant Branch (AGB) progenitor's slow wind (e.g., 
Frank & Mellema 1994). A shortcoming of the hydrody- 
namic modeling is that often several sets of parameters 
based on differing geometries are capable of reproducing 
the same morphology. Observations of spatially resolved 
kinematics of a PN can help differentiate among models. 

We have been studying the physical structure of PNe 
using both imaging and high-dispersion spectroscopic ob- 
servations, in order to understand their formation and 
evolution. Thus far, we have analyzed the Owl Neb- 
ula (NGC 3587) in Paper I (Guerrero et al. 2003) and 
NGC 7662 in Paper II (Guerrero et al. 2004). The Owl 
Nebula is a triple-shell PN; its "owl face" morphology is 
caused by bipolar cavities excavated by a fast wind that 
ceased some time ago, and its kinematics indicate that the 
nebular material is backfilling these cavities. NGC 7662 
is also a triple-shell PN; its bright main nebula consists 
of two elliptical shells with simple structures in Ha and 
[O III] , although the outer rim of the main nebula is in- 
terspersed with [N II]-bright knots. These two PNe were 
chosen for our study because their morphologies have been 
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than the optimal —88 °C at which it normally operates). 
This higher operating temperature produced a substan- 
tially greater number of warm pixels and a higher dark 
current. To process the observations we first flagged and 
interpolated the warm pixels using the hot pixel lists avail- 
able at STScI 3 , and then proceeded with the standard 
routines in the STSDAS package within IRAF. These im- 
ages are not shown in this paper because the F658N im- 
age did not detect any nebulosity and the F656N image 
does not show any small-scale features unresolved in our 
ground-based Ha image in Figure 1. The F656N image 
has only been used to flux-calibrate our ground-based Ha 
image. The total Ha flux from NGC 1360 is -3.1 xlCT 10 
ergs cm -2 s _1 . 

2.2. High- Dispersion Spectroscopic Observations 

High-dispersion spectroscopic observations of NGC 1360 
were made using the echelle spectrograph on the 4 m tele- 
scope at Cerro Tololo Inter- American Observatory (CTIO) 
on 2001 December 25. The 79 line mm -1 echelle grat- 
ing and the long-focus red camera were used to achieve a 
reciprocal dispersion of 3.4 A mm^ 1 . The spectrograph 
was used in the single-order, long-slit mode to allow the 
coverage of the Ha and [N n] AA6548,6584 lines along a 
slit length of ~3'. The data were recorded with the SITe 
2K No. 6 CCD with a pixel size of 24 /zm, which corre- 
sponds to ft' 254 pixel -1 along the spatial direction and 
3.7 km s _1 pixel" 1 along the dispersion direction. The 
slit width was l'.'O, and the resultant instrumental FWHM 
was 12 km s _1 at the Ha line. 

The echelle observations were made along the major axis 
of NGC 1360 at the position angle (PA) of 27°. Three slit 
positions were used to cover the full extent of the major 
axis, as marked in Figure 1. The integration time at each 
slit position was 600 s. The angular resolution, determined 
by the seeing, was — 2'.'0. A mosaic of the Ha echellograms 
is shown in Figure 2, alongside an Ha image oriented such 
that the major axis is parallel to the spatial direction of 
the echellogram. 

3. RESULTS 

3.1. Morphology 

As shown in the Ha image (Fig. 1), NGC 1360 ap- 
pears elliptical with its major axis oriented along PA — 
30°. It lacks a sharp outer boundary because the sur- 
face brightness declines with distance until it drops be- 
low the background noise. In order to estimate the size 
of NGC 1360, we determined the background noise in the 
Ha image and defined the nebular edge to be at 3a above 
the background. Using the threshold of 3er = 5.3 x 10~ 16 
ergs cm~ 2 s _1 arcscc~ 2 , NGC 1360 extends 9' x 7', as 
shown by the contour in Figure 1. At a distance of 
350j^5o°° pc as determined from the Hipparcos measure- 
ments (Acker et al. 1998), the linear size of NGC 1360 
would be 0.9 x 0.7 pc. 

While the overall surface brightness distribution of 
NGC 1360 falls off radially, there are depressions of —20% 
within 30" of the central star and in the eastern hemi- 
sphere. The spatial distribution of these surface bright- 
ness depressions do not appear to be random. The central 

3 http://www.stsci.edu/instrumcnts/wfpc2/wfpc2_hotpix.html 



depression extends into the eastern hemisphere and curves 
southward, while another distinct depression curves from 
the east to the north of the region surrounding the central 
star. The cause of these depressions is either lower gas 
densities or a higher extinction. If the latter is true, the 
non-random distribution of the depressions suggests that 
the obscuring material is internal to the nebula. A quan- 
titative analysis of the surface brightness of NGC 1360 is 
presented in §4. 

At the northern end of the major axis there exist sev- 
eral nebular knots extending from 240" out to 350" from 
the central star. Near the southern end of the major axis, 
there arc also knots (marked by an arrow in Fig. 1), but 
they are fainter and detected at just the 3u level. The faint 
southern knots, the central star, and the bright northern 
knots are not aligned along a straight line; instead, they 
form a slightly curved arc with the star at the apex. Such 
a configuration has been observed in other PNe and sug- 
gested nebular motion through the interstellar medium. 

3.2. Kinematics 

The echelle observations at the three consecutive slit 
positions have been mosaicked into an echellogram span- 
ning the full major axis of NGC 1360 and are presented in 
Figure 2. This Ha echellogram shows a position- velocity 
ellipse, indicating an expanding shell in the nebula. Fur- 
thermore, the position-velocity ellipse is tilted, with the 
southwestern and northeastern tips blue- and red-shifted, 
respectively. For an expanding ellipsoidal shell, the tilt of 
the position-velocity ellipse results from the inclination of 
the shell's polar axis with respect to the sky plane. As 
the position-velocity tilt occurs along the major axis of 
NGC 1360, its ellipsoidal shell must be prolate in shape, 
with its southwestern tip oriented toward us and its north- 
eastern tip away from us. 

The velocity profile of the Ha line near the central star 
shows two broad components (see Fig. 3). The approach- 
ing component has a heliocentric velocity of Vhci = 25 
km s _1 and a width of FWHM = 30 km s _1 , while the 
receding component has Vhci = 73 km s _1 and FWHM = 
37 km s _1 . These widths are much larger than the thermal 
width of the Ha line at 10 4 K, FWHM = 21.4 km s"\ indi- 
cating the existence of significant velocity variations along 
the line of sight within the nebular shell. 

The average velocity of the approaching and receding 
sides of the nebula, Vhci = 49 km s _1 , is adopted as the 
systemic velocity of NGC 1360. This result is in good 
agreement with the previous measurement of Vhci = +49.1 
km s" 1 by Bohuski & Smith (1974) using a Fabry-Perot 
scanner for a 60" aperture near the center of NGC 1360. 
The line-splitting at the center of a PN is usually adopted 
as twice its expansion velocity so that the apparent expan- 
sion velocity of NGC 1360 would be 24 km s -1 ; however, 
as we discuss in §4.1, this is an oversimplification of the 
expansion of NGC 1360. 

In the Ha line, both the nebula and the knots are de- 
tected; however, in the [N II] line, only the knots are de- 
tected. This is illustrated by the Ha and [N II] lines of 
the northeastern end of NGC 1360's major axis shown in 
Figure 4. The low [N II] /Ha ratio of the nebula is ex- 
pected because the central star of NGC 1360 has an cf- 
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fcctivc temperature of 110,000 K (Napiwotzki 1999) and 
the ionization stage of N in the nebula is higher than N + . 
The high [N II] /Ha ratios of the knots suggest that they 
may be FLIERS, i.e., fast, low-ionization emission regions 
(Balick ct al. 1993). Figure 4 shows that the innermost 
knots detected within the slit are at ~260" from the central 
star and have a radial velocity of Vhel ^100 km s _1 corre- 
sponding to a systemic velocity of ^+50 km s _1 . The ra- 
dial velocity of the knots increases outward by 0.5 km s _1 
arcsec -1 . These knots clearly move faster than the shell 
expansion, confirming their identification as FLIERs. 

4. DISCUSSION 
4.1. Spatio- Kinematic Structure of NGC 1360 

The Ha image and echelle spectra of NGC 1360 can be 
used to construct spatio-kinematic models. In contrast 
to the identifiable thin shells in the Owl Nebula (Paper 
I) and NGC 7662 (Paper II) that can be modeled easily, 
the shell of NGC 1360 has both large velocity widths and 
irregular brightness variations in the echelle line image, 
which require simplifying assumptions to make the mod- 
eling tractable. We have assumed that NGC 1360 consists 
of a series of concentric thin ellipsoidal shells and the shell 
expansion velocity increases outwards. Consequently, the 
envelope of the position- velocity ellipse (marked in Fig. 2a) 
corresponds to an outer thin ellipsoidal shell that may be 
identified with the elliptical isophotes with a matching spa- 
tial extent in the direct image (marked in Fig. 2b). This 
outer thin shell can then be modeled using the methodol- 
ogy outlined by Guerrero et al. (2003), by assuming that 
the shell expands homologously (i.e., a Hubble law). 

We model this shell of NGC 1360 as a prolate ellipsoidal 
shell with its polar axis tilted with respect to the line-of- 
sight by an inclination angle i. Its polar and equatorial 
radii are r p and r e , and its polar and equatorial expansion 
velocities are V p and V e , respectively. For the assumed ho- 
mologous expansion, r p /r c is equal to V p /V c , and the kine- 
matic age is simply r e /V e or r p /V p . A grid of models are 
constructed for a range of inclination angles i, axial ratios 
r p /r e , and equatorial expansion velocities V e . These mod- 
els are then used to synthesize nebular isophotes and to 
generate position-velocity ellipses for visual comparisons 
with the observations. Our best-fit model is marked in 
Figure 2, and has the following parameters: i = 60°, V c = 
35 km s" 1 at r c = 130", V p = 70 km s" 1 at r p = 260", 
and ikin = 10,000 yr. Note that the inclination angle and 
the axial ratio are correlated so that lower values of i yield 
higher ratios of V p /V e . Despite this fact, the inclination 
angle and the axial ratio are both well constrained. The 
former is within ±10° from the optimal value, while the 
latter is within ±0.25 from V p /V e = 2.0. The expansion 
velocities and kinematic age, on the other hand, are not 
as well constrained. For i = 60°, the uncertainty in V e is 
±10 km s _1 , corresponding to about ±30%. 

Bear in mind that this model describes only an outer 
layer of the nebular shell of NGC 1360. The bulk of the 
nebular material is interior to this layer and expands more 
slowly. This densest shell material is responsible for the 
peaks in the velocity profile. For an inclination angle of 
60° , the velocity split near the nebular center roughly cor- 
responds to twice the equatorial expansion velocity of the 
densest nebular material in NGC 1360, 24 km s _1 (see 



Fig. 3 and §3.2). However, the detailed geometry and 
kinematics of this dense material cannot be more precisely 
determined because the peaks of the line profiles do not 
form a well-defined position-velocity ellipse in the echelle 
image. 

The FLIERs in NGC 1360 are projected along the ap- 
parent major axis of the nebula; thus, it is reasonable to 
assume that the FLIERs are distributed along the polar 
axis with an inclination angle of 60°. The apparent ex- 
pansion velocities of the FLIERs, >50 km s , then corre- 
spond to deprojected expansion velocities of >100 km s _1 . 
These expansion velocities are typical of those observed in 
FLIERs (Balick, Preston, & Icke 1987). 

4.2. Physical Structure of NGC 1360 

We can use the spatio-kinematic model of NGC 1360 de- 
duced in §4.1 along with the observed Ha surface bright- 
ness distribution to study the density structure within 
this nebula. The observed Ha surface brightness profiles 
along several PAs are plotted in Figure 5. The depres- 
sion in surface brightness seen in the eastern hemisphere 
of NGC 1360 is clearly visible in the profiles along PAs be- 
tween 105° and 165°. In contrast, the surface brightness of 
the western hemisphere of NGC 1360 shows smoother pro- 
files, without obvious "bumps" or "dents." We have carried 
out the analysis of the density structure of NGC 1360 us- 
ing the sections of the nebula between PA 225° and 345° 
because their uniformity best represents the overall distri- 
bution of material in the nebular shell. 

The profiles of the western regions of NGC 1360 in Fig- 
ure 5 exhibit a Gaussian shaped brightness distribution 
peaking at 40"-60" and declining smoothly with increasing 
radius. The smooth shapes of the surface brightness pro- 
files indicate that the nebula does not contain any regions 
with sharp density or extinction variations in the west- 
ern hemisphere. These surface brightness distributions, 
in conjunction with the geometric information determined 
from the spatio-kinematic modeling in §3.2, can be used to 
derive the nebular density distribution. We have modeled 
the surface brightness profile assuming a Gaussian density 
distribution and fitted this model to the observed Ha sur- 
face brightness profiles between PA 225° and 345°. The 
best-fit model peaks at 35" with a FWHM of 60" along 
the equatorial plane and peaks at 70" with a FWHM of 
100" along the polar direction. The best-fit models are 
overplotted in Figure 5: good agreements are obtained 
from PA 225° to 345°. Using the flux calibrated surface 
brightness we further determined that the peak density is 
~130 H-atom cm -3 , and that the density at r c = 130" 
is ~0.1 H-atom cm" 3 . The low density at the outer edge 
of the nebula implies that it is beginning to merge into 
the interstellar medium. The total mass in the model, 
~0.4 M©, is estimated assuming that the surface bright- 
ness depressions in the eastern hemisphere are caused by 
internal dust absorption. If the lower surface brightness is 
caused by a lack of material in the shell, the mass would 
be over-estimated, but this error is small as compared to 
the uncertainty in mass due to the large error bar in the 
distance to NGC 1360. 

5. SUMMARY AND CONCLUSIONS 
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We have obtained images and high-dispersion long-slit 
spectra of NGC 1360 in order to study its spatio-kinematic 
structure. We have modeled an outer layer of the shell as 
a prolate ellipsoid and find a major-to-minor axial ratio 
of 2:1 and a polar axis tilted by 60° with respect to the 
line of sight. The kinematic age derived from this model 
is -10,000 yr. 

We have found that the observed morphology and sur- 
face brightness profiles of NGC 1360 can be modeled well 
by a thick shell with a Gaussian radial density profile. 
This smooth density distribution implies that there is no 
sharp inner edge as expected from compression by a fast 
stellar wind. If a fast wind was present in the past and 
compressed the AGB wind into a thin shell, that shell 
has since diffused. The low nebular density, <130 H-atom 



cm~ 3 , suggests that NGC 1360 is an evolved PN; its low- 
density (0.1 H-atom cm~ 3 ) envelope has begun to merge 
with the interstellar medium. 

The FLIERs observed in NGC 1360 are intriguing. The 
kinematic age of the FLIERs, based on their de-projected 
distances from the central star, >0.5 pc, and their de- 
projected expansion velocities, > 100 km s _1 , is ~5,000 
yr. This kinematic age is about 1/2 that of the PN itself. 
NGC 1360 is one of the very few cases of evolved PNe in 
which collimated outflows are observed (e.g., Corradi et 
al. 1999). 

We thank Jin-chung Hsu and Stefano Casertano for their 
help in addressing the warm pixel problem of the HST 
WFPC2 image. 
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Fig. 1. — An Ha image of NGC 1360 displayed with two different intensity contrasts to highlight the faint outer regions (left) and the 
structure in the central regions (right). The contour in the left panel corresponds to the 3(7 level above the background. The slit positions 
used in the echelle spectroscopic observations are marked on the image. The faint southern FLIERS are marked with arrows in both panels. 



Fig. 2. — (left) Ha echellogram along the major axis of NGC 1360. The gcocoronal Ha line is marked, (right) Ha image of NGC 1360 
oriented to match its major axis with the echelle slit position for direct comparison. The jagged contours overlaid on the echellogram and 
direct image have matching spatial extent and correspond to isophotes detected at 4c and 7a levels, respectively. The smooth ellipses overlaid 
on the echellogram and direct image are based on our best-fit spatio-kincmatic model. The bright FLIERs northeast to the nebula are visible 
in both panels. 



Fig. 3. — Velocity profile of the Ha line within the central region of NGC 1360. The position of the geocoronal Ha line is marked. The 
systemic velocity adopted for NGC 1360 is also marked. 



Fig. 4. — An expanded view of the Ha and [N o] echcllograms at the northeast end of NGC 1360 to show the FLIERs. The spatial position 
is measured relative to the central star. The Ha line shows both the nebular shell and the FLIERs, while the [N II] line shows only the 
FLIERs. The FLIERs exhibit a velocity gradient along the slit visible in both lines. 



Fig. 5. — Surface brightness profiles (solid lines) of the Ha line extracted at 30° intervals. The peak at the zero position corresponds to the 
central star, and the position angle is marked for the profile on either side of the star. The dashed curves are the profiles from our best-fit 
model. 
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